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POGO  STABILITY.  RELIABILITY  AND  PARAMETERS  ANALYSIS 

Wang,  Qisheng:  Gau.  Wanyong:  Gu.  Yongchun;  Zhang.  Jitony; 

and  Li .  Xianshan 

Abstract  POGO  instability  is  one  of  the  most  important  dynamical 
problems  of  space  vehicles.  The  mechanism  of  the  POGO  instability  is 
briefly  clarified  by  a  simple  free-free  system  which  consists  of  two 
masses,  a  spring,  and  a  damper,  followed  by  POGO  stability  matrix 
algorithm.  POGO  stability  single-transfer  method,  and  POGO  stability 
estimation  method-  The  adaptability  of  these  POGO  stability  analysis 
is  illustrated  through  an  example.  After  important  parameters 
related  to  POGO  instability  are  discussed,  the  reliability  analysis 
of  POGO  stability  and  its  criteria  are  also  introduced. 

(This  paper  was  received  on  September  6.  1984.) 

I.  Mechanism  of  POGO  Instability 

The  phenomenon  of  POGO  not  only  has  adverse  effect  on  the 
passenger-carrying  space  vehicles,  it  is  also  an  important  element 
in  structural  design.  The  high  pressure  vibration  in  the  propellent 
line  caused  by  POGO  can  result  in  the  degradation  of  the  performance 
of  the  propulsion  system,  and  sometimes  even  pre-mature  engine  shut¬ 
down  [1.2] . 

In  the  70 's.  the  suppression  of  POGO  instability  was  viewed  as 
one  of  the  most  important  criteria  in  the  design  of  gas-1 iquid 
structural  system  and  manned  space  vehicle.  In  the  80 ‘s.  it  was 


I 


still  a  part  of  the  research  of  space-traveling  related  structural 
,  dynamics  (31. 

Similar  to  flutter  and  galloping  in  the  gas  elasticity,  POGO  is 
an  important  problem  in  the  fluid-solid  coupled  dynamics.  Flutter 
anu  galloping  are  the  dynamic  problems  caused  by  the  flowing  of 
fluid  exterior  of  a  solid,  while  POGO  is  the  dynamic  problem ’caiiSed 
by  the  flowing  of  fluid  interior  of  a  solid.  Since  the  structural 
vibration  is  a  three-dimensional  problem,  and  the  pipeline  system  in 
real  space  consists  of  pipelines  in  all  three  dimensions  which  is 
further  controlled  by  the  propulsion  of  the  engine  and  other 
controlling  forces.  POGO  is  the  dynamic  problem  of  a  large-loop 
coupling  system. 

POGO  is  the  study  of  structural  instability  of  the  rocket  body 
aitructural  system  and  the  longitudinal  coupling  of  the  propulsion 
system.  Under  the  condition  of  POGO,  longitudinal  self-stimulating 
non-convergent  vibration  of  the  space  vehicle  is  caused  which  is 
similar  to  the  ever- jumping  "pogo  stick’’  and  hence  was  given  the 
name. 

The  stability  of  POGO  was  already  discussed  in  literatures  (4] 
to  [7] .  In  the  past,  the  mechanism  of  POGO  was  mostly  discussed  in 
words  or  through  complicated  calculation.  Similar  to  the  mechanism 
of  flutter  [8] ,  POGO  can  be  explained  by  the  system  of  two  free-free 
masses(Mi,  M2)  -  damper(C2)  -  spring(K2).  figure  1.1b,  we  have 


M^X^+CtU^-X,HKt{Xt-X,)=‘0 


let  d^X,-X„  oi=/C.(^,+A/.)/(M,M.) 

X 


(1.2) 


(1.3) 


2£,o),-C.(M,  +  Ai.VCM.A/.) 

Assuming  that  the  pulse  thrust  FjL(t)  is  proportional  to  the  relative 
deformation  velocity'^,  then 

r=F.(/)/.4/,=2^<K».^  4) 

where  Ce  is  related  to  many  factors  of  the  propulsion  system.  Using 
the  above  equation,  we  have 


if  the  initial  perturbation  is  3„  then  the  solution  can  be  , 

and  we  have 


S*  +  20»Un-Cc)S  +  Ql^O 


the  roots  are 


(1,6) 


6’i,i  = -(£-."fc)<a-±/o),\/ l-(  .  ^ 

and  the  solution  is  therefore. 

<3=<J.EXF[-(S.-i.)o.Ocos(o)W  +  0.) 
where  ' 

0,  is  a  constant . 

If  equation  (1.8)  is  plugged  into  equations  (1.1)  to  (1.3). 
then  the  reverse  motion  of  X|  and  '/.o  can  be  realized.  If  the  actual 
structural  modular  damping  ratio  f*. i.s  a  positive  dampina,  then  when 

y 


then  the  unstable  divergent  vibration  due  to  the  negative 
damping  of  the  system  will  occur.  In  this  case,  the  system  will 
absorb  energy  and  the  source  of  energy  is  the  work  of  pump  and 
engine.  This  is  similar  to  the  flutter  of  wings  of  aircraft,  even 
though  the  engine  provide  the  flight  speed  to  the  aircraft,  not  all 
the  flight  speed  will  cause  unstable  flutter  of  the  wings.  For  POGO 
systems,  only  under  certain  coupling  conditions  v/hich  render  large  f,: 
can  unstable  pogo  occur  and  the  extent  of  instability  if  determined 
by  the  magnitude  of  £ .  The  magnitude  of  £«  is  further  determined  by 
structure,  pump,  and  the  closeness  of  the  gain  of  the  engine  and  the 
characteristic  frequency  of  the  propellant  line.  The  latter  is  an 
important  design  criterion  of  POGO  (see  section  V.3).  If  coupling  is 
removed  by  some  measures.  £».can  be  reduced  and  system  POGO  can  be 
avoided. 

As  mentioned  above.  POGO  of  space  vehicles  is  a  problem  of  the 
large-loop  coupling  dynamics.  Since  the  structure  of  space  vehicles 
consisted  of  many  modules,  many  propellants,  many  pipelines,  and  the 
simultaneous  working  of  pump  and  engine,  some  local  vibration  can 
also  couse  POGO.  theref  re.  the  dynamic  equations  for  many  parts  can 
be  generated  (see  section  II).  Then  the  dynamic  polynomials  of  the 
entire  POGO  system  can  be  tabulated  and  then  the  calculation  of  the 
POGO  stability  of  the  entire  system  can  be  carried  out  (see  section 
III) .  However,  since  the  physical  significance  is  hard  to  detect  in 
the  initial  stage  of  calculation,  the  POGO  single-transfer  and 
estimation  methods  should  be  used  (see  section  IV)  to  estimate  the 
POGO  stability  of  the  vehicle  in  a  more  efficient  way.  Nevertheless, 
the  relative  error  should  be  noted  in  the  estimation  procedure. 


There  are  many  factors  affectina  POGO.,  If  the  major  factors  ar 
identified,  the  requirements  for  the  system  design,  test,  planning, 
and  quality  control  can  be  outlined  (see  section  V).  Furthermore, 
there  are  some  uncertain  parameters  in  the  POGO  stability  analysis 
and  the  reliability  analysis  (section  VI)  and  the  POGO  reliability 
requirement  which  suppresses  POGO  (section  VII)  should  be  carried 
out . 

II.  The  Related  Dynamic  Equations  of  a  POGO  System 

The  typical  block  diagram  of  the  simplified  coupling  of 
propulsion-structure  system  is  shown  in  figui'e  4.1.  Four  major  sub¬ 
systems  are  shown;  structural  svstem  of  the  vehicle,  propellant 
pipeline  system,  pump  system,  and  the  thrust  c.hamber  system.  The 
description  of  these  systems  is  followed. 

(I)  Dynamic  equations  of  fne  vehicle  structural  system 
The  dynamic  equations  for  tne  ith  degree  modular  state  are 


x,(e)-c.,(srrM 

)  -  .SO-,((.>)/(6  =  --ZC.co.S  -  co} ) 


1  _ 


(2.1.1) 

(2.1.2) 

(2.1.3) 

(2.1.4) 


where  iV/(,  &),»  0((;)»  and  Gi(c)  are  the  modular  damping  ratio, 

macroscopic  mass,  circular  frequency  of  the  ith  modular  state,  the 
modular  displacement  of  the  jth  point,  general  modular  displacement 
of  the  engine,  and  structural  gain.  Ffji  is  the  external  force 
applied  on  point  j  of  the  roc.kei  ix-dy.  The  physical  displacement  of 
the  jth  point  on  the  structure  is 


(o 


^(/)-  Z  =  S  (2.1. 

<»l  <st 

Where  is  the  macroscopic  coordinate:  The  form  of  the  typlcsl 

external  force  F{j)  is  explained  below. 

1.  The  pulse  thrust  F(l)  created  by  the  pulse  pressure  Pc  of 
the  combustion  chamber  on  the  engine  is 


Fiiy-=A,kC/F.--^S,l\  (2.1.6) 

where  Se^A^^^f.  and  A^h  and  Cf  are  the  cross  section  of  the  neck  of 
the  combustion  chamber  and  thrust  coefficient. 

2.  The  forces  F(2)  and  F(3)  applied  on  the  structure  as  a 
result  of  the  pulse  pressure  due  to  elbow  or  pipelines  before.'^af tei 
a  pump  and  momentum  change  are  (see  figure  2.3) 


“  AuPu  —  pQu(^2Qi// At) -i- JCi)  (2.1.7) 

F(  3  )=  -AoPo  -  pQoi2QolAa  -  F*) 

(2.1.8) 

where  Aij,  Py,  Qu^.  Qy  are  the  inlet  cross  section,  pulse  pressure, 
steady  state  and  pulsed  volume  flow  of  the  upstream  flov?:  and  Ar, . 
PD'  Qot  Qd  corresponding  quantities  of  the  downstream  flow. 

Xi  and  fi  are  the  transverse  and  longitudinal  velocity  at  the  elbow 
or  pump. 

3.  Due  to  change  in  cross  sectional  area,  the  change  in 
pressure  and  momentum  result  in  an  applied  force  on  the  structure 
F(4)  (see  figure  2.4)  is 


F(4)- -AuPu  :  AnPo-2pQ( 


(2.1.9) 


4.  The  applied  force  F(5)  as  a  result  of  the  viscosity  of  the 


fluid  in  a  straight  line  (see  figure  2.1)  is 


-  For  SLAX,) 


(2.1.10) 


where  A.  R.  Z.  L.  Jt,  are  the  corresponding  cross  sectional  area, 
resistance,  impedance,  inertia,  and  velocity. 

5.  The  applied  force  created  iiy  the  pulsed  flow  Qp  of  the 
propellent  on  the  bottom  of  the  fuel  tank  on  the  center-of-mass  of 
the  fuel  tank  with  respect  to  the  structure  is 


F(6)=pA,^a.  (2.1.11) 

where  A,;  is  the  height  of  the  liquid  in  the  fuel  tank. 

6.  The  applied  force  F(7)  created  by  the  pulsed  pressure  .^t  at 
the  opening  (tb)  at  the  bottom  of  the  fuel  tank  on  the  bottom 
structure  is  (6] 


F(7)=^.P, 


(2.1.12) 


where 


F,  —  ^^ph,ii>Kin^K~  (''.1.13) 

Qi,  =  Qir  AiXitb)~Q,-\rAx'^(^K(Jb)iK  (2.1.14) 

A’»i 

where  4>kU)>  and  Lt  are  the  modular  displacement  of  the  center- 

of-mass  of  the  liquid  reservoir,  the  modular  displacement  of  the 


bottom  of  the  tank,  arid  inertia  of  the  liquid  in  the  tank.  A^  is  the 
area  of  the  opening,  and  Qr  and  jt  are  the  relative  and  true  flow 


7.  The  force  F(8)  acting  on  the  structure  of  the  rocket  body  as 


a  result  of  the  local,  longitudinal  motion  of  the  pump  is 


(2.1.15) 

where  Jif(tp)  is  the  acceleration  of  the  pump  component  of  the  rocket 
structure.  Xp  is  the  acceleration  due  to  local  longitudinal  motion 
of  the  pump,  and  Mp  is  the  mass  of  the  pump  component. 

(II)  Dynamic  equations  of  the  pipeline  sections 
1.  straight  pipeline  sections  (see  figure  2.1) 

Assuming  that  the  steady  state  velocity  in  the  pipeline  is  much 
less  than  the  speed-of -sound,  then  the  effect  of  Mach  number  can  be 
ignored  and  the  following  relationship  holds 


chO 

-jZshO  ' 

'■^ARshO  '| 

jOshO 

chO  j 

X, 


(2.2.1) 


where  fl'=6''r'Z/(5L)j  r=ljar  :  1  is  the  length  of  the  pipeline. 

L=pllAt  Z=R-¥SL  :  a^  IS  the  equivalent  speed-of -sound  '.see 
discussions  in  section  Vfo)). 

If  the  effect  of  viscosity  can  be  ignored  than  R  is 
approximately  zero  and  the  above  equation  has  the  following  form 


cA(t«) 

£.«A(t«)'| 

r 

io.j 

Tshivs') 

<  ~r~ 

cA(«) 

< 

(2.2.2) 


If  the  effect  of  compressibility  can  be  ignored,  then 
approaches  infinity  and  the  above  equation  is  further  simplified  as 


Wo  1X1 


(2.2.3) 


2.  If  the  pipeline  system  is  complicated  and  the  gas  content  of 
the  liquid  and  the  gas  bubble  generated  during  flow  of  the  liquid 
can  not  be  estimated  easily,  then  the  calculation  of  equivalent  •  : 
speed-of-sound  is  made  very  difficult  and  the  modular  parameters  of 
the  pipeline  system  can  be  estimated  by  the  trial  method.  For 
example,  when  the  two  ends  of  the  pipeline  connecting  the  bottom  of 
the  liquid  reservoir  and  the  pump  is  closed/opened,  then  we  have 

(S*  +  2i,a>,S+a))q,  =  A,Pu/M, 

(2.2.4) 

When  the  two  ends  are  closed/closed,  then 

(2.2.5) 

where  Cu  o/.  Mu  qi-,  •  are  the  modular  damping  ratio,  circular 
frequency,  macroscopic  mass,  and  macroscopic  coordinate  of  the 
pipeline . 

3.  The  equation  for  the  pressure  retainer  (or  connecting 
between  branches)  is  (see  figure  2.2) 

iQo>  L-y,  (2.2.6) 

where 

ya  =  SCt/(.LaCaS'  +  R,C,S  + l)t  ^2  2  7) 

Lt  =  plt/A,f  +  A»  =  nn.rl  (2.2.8) 

r^tPiKyP’)  (isothermal  process) 

C.=y/ivP)  =  \ 


d 


(2.2.9) 


‘  liyP)  (adiabatic  process) 

R,  =  ^i^PjKi  K=-2pgKCiA,y  (2.2.10) 

where  and  r^  are  the  length  and  radius  of  the  manifold  of  the 
pressure  retainer,  n^  is  the  number  of  holes,  v  is  the  ratio  of  the 
specific  heat  of  the  gases  in  the  pressure  retainer,  p is  the  mass 
density  of  the  liquid.  is  the  outlet  coefficient  of  the  short 

tube.  Pq  and  Vq  are  the  initial  pressure  and  volume  of  the  pressure 
retainer.  P  is  the  meta-stable  inlet  pressure  of  the  pressure 
retainer  during  the  flight.  AP  is  the  pulse  pressure  difference  of 
the  ends  of  a  short  tube  when  resistance  is  measured,  and  g  is 
the  gravitational  acceleration. 

(Ill)  Dynamic  equations  of  a  pump  system  (see  figure  2.3) 


(Pd  )  SC pZ p  "Zp~\tJ*(/)  r~^liiZp  ~ S Lp.-loif  Jl p ) 

io„rL  -scp  1  JiQ„rL  Jii'pi 


(2.3.1) 


where  (m+1),  Cp.  Lp,  Rp,  Zp=Rp+3Lp  are  the  dynamic  gain  of  tlie  pump. 

softness  of  the  gas  clock,  inertia,  resistance  and  impedance.  When 

•  • 

(m+l)“i  and  =  the  above  equation  has  the  following  form 


(Pd  )_rl-rSCpZp  -Zp~i(Pu  )  (2.3.2) 

(Q«3  L  -SCp  1  Jc?,,) 

when  Zp-0  (only  for  the  composite  parameter  Cp) ,  the  above  equation 
is  transformed  to 


0 

1 


f  ^ 


(2.3.3) 


(IV)  Related  equations  for  the  engine  (see  figure  2,5) 

1.  equation  for  the  after-pump  pipeline  sections  (not 
considering  the  compressibility  of  the  liquid) 


ri 

(2.4. i) 

iQj 

Lo  1  Jio„) 

where  2d-Rd+SL£).  and  Rj).  Lj)  are  the  resistance  and  inertia  of  the 
after  pump  pipeline  section,  and  i'Q  is  the  pipe  wall  velocity. 

2.  equation  for  the  nozzle 

When  the  compressibility  of  the  liquid  in  the  nozzle  section  i 
not  considered  (similar  the  treatment  above),  we  have 


Pj-Pc+ZjQj-AjRjfi  (2.4.2) 

where  Zj-Rj+SLj  and  A j .  Rj,  L j .  and  are  the  cross  sectional  area 
resistance,  inertia,  and  velocity  of  the  nozzle  section. 

3.  combustion  chamber 

From  the  combustion  thermodynamics,  the  relation  between 
combustion  chamber  pulse  pressure  and  two  pulse  propellant  flow, 
rates  Oj  and  is 


Pc  =  Zc(0/  +  0*) 


(2.4.3) 


where  (1  +  rcS),  Hg-C*/ ( Ath9)  •  ^th*  ^c 


neck  cross  sectional  area,  characteristic  velocity,  time  delav.  and 


impedance. 


III.  POGO  Stability  Matrix  Alaorithm 


There  are  several  methods  which  consider  the  stability  of  the 


system.  Most  of  these  methods  require  the  solution  of  the  close- 

loop/open-loop  transfer  coefficients  and  the  zeros,  polar  points  and 
the  numbers,  and  then  the  stability  of  the  entire  close-loop  can  be 
analysed. 

The  transfer  coefficient  for  some  systems  can  be  deduced  from 
the  trial  method.  However,  it  is  very  difficult  to  measure  the 
transfer  coefficient  of  a  close-loop  svstem  from  an  unstable  c>  se- 
loop  test.  Similarly,  there  are  stable  and  unstable  open-loop 
systems  and  the  determination  of  their  transfer  coefficients  is 
equally  difficult.  Therefore,  it  is  very  important  to  evaluate  the 
system  characteristics  (close  or  open  loop)  based  on  the  composite 
calculation  of  the  refined  model  of  many  sub-systems. 

If  a  pseudo  element  or  pseudo  block  is  introduced  into  the 
series  of  equations  of  a  close-ioop  system  or  its  block  diagram,  the 
characteristic  equation  of  a  c l';'Se-looD  or  the  transfer  coefficient 
and  the  zero  polar  point  of  an  ioen-loop  system  can  be  determined  in 
a  unified  and  convenient  way.  Df^.ause  the  location  of  the  opening  is 
not  unique,  the  introduction  -jf  a  pseudo  block  is  not  unique  either. 
The  selection  of  the  location  for  opening  or  testing,  therefore, 
should  be  based  on  experience.  Of  course,  it  is  best  if  the  pseudo 
block  can  be  placed  on  the  main  loop  of  a  close-loop  system. 

Assuming  that  there  is  no  e.xterior  disturbance,  we  can  obtain  a 
series  of  equations  for  the  n  variables  of  the  n  equations  of  the 
POGC  system  as  explained  in  section  II.  To  illustrate  this  point, 
assume  (or  deduce  from  simplification)  that  n=3.  the  followina 


series  of  equations  .( in  matrix  form)  correspondinq  to  block  diaaram 


3 . 1  is 


*  0 

Dc(S) 

-McdS)  - 

rr 

DniS) 

-MniS) 

0 

.  -K„ 

0 

1  . 

r 

■LAiSmX}-^0 


(  3 , 0  .  i  ) 


where 


*  0  Dc(5) 

-McCsy 

rr 

C^(5)]  = 

DniS)  -MniS) 

0 

{A-}=. 

.  -Ko  0 

1*  , 

T' 

S  is  the  LaPlace  variable.  M(S)  and  D(S)  are  the  corresponding 
polynomials  of  the  numerator  and  denominator. 

Even  though  there  are  many  analytical  methods  for  the  stability 
analysis,  the  basic  nature  of  these  methods  is  similar:  namely,  the 
system  is  stable  if  the  real  part  of  the  root  of  the  close-loop 
characteristic  equation  if  negative,  otherwise  the  system  is 
unstable.  Two  analytical  methods  emanate  from  this  point:  the  close- 
loop  and  open- loop  analytical  method. 

1.  Direct  root-finding  from  the  characterist ic  equation  of  the 
close- loop*  system 

From  the  real  and  imaginary  parts  of  the  root  of  the 
characteristic  equation  !A(S)!-0.  one  can  determine  the  stability  of 
a  POGO  system  in  any  frequency  range.  Under  the  perturbation-free 
condition,  the  determinant  of  the  variables  of  the  nth  order 
polynomial  is  zero  which  is  the  characteristic  system  equation.  If 
there  is  a  pseudo  block,  let  K]-)=i  and  one  would  have 


M( 5)  I  -  Dc( S)  D/,(  S)  -  S)M„( S)  =  0 


(3.1.1) 


The  form  of  root  of  the  characteristic  equation  is 


where 


Sk~Si,  ±  jSki—o)i,U±ja>i,y/l—Cl 


(3.1.2) 


(3.1.3) 


where  and  Sj^j  are  the  real  and  imaginary  parts  of  the 
characteristic  root:  wj,.  and  Ck  are  the  typical  frequency  and  damping 
ratio  of  the  POGO  system  in  question. 

If  the  damping  ratio  of  the  nth  modular  state  of  the  original 
real  structure  is  and  if  the  damping  ratio  of  certain  modular 
state  is  changed  from  S,-  to  Sc-  (while  other  damping  ratios  are 

either  kept  constant  or  changed  with  the  same  ratio),  then  v/e  should 

» 

have  at  least  one  neutral  stable  root  (if  Sp^j^-O)  and  the  width  A 
(dB)  of  the  stable  gain  of  the  system  is 


8)=^ 20 LogMJCc).  (fc>0)  (3  14) 

If  Ce<Q.  this  would  indicate  that  the  system  is  already  a  stable 
device.  The  system  is  stable  if  A -O.  whereas  the  system  is  un.stable 
if  A<;-0. 

2.  Open- loop  analytical  method 

Taking  advantage  of  the  pseudo  block  method  and  let  Kj;,-0.  one 
can  obtain  the  characteristic  equation-polar  point  equation  of  the 
open- loop  system 


5 )  =  M( 5 ,  =  0)  I -=  DcC S )D„( 5 )  =  0 


(3.2.1) 


the  magnitude  and  number  of  the  roots  of  the  open-loop 

characteristic  equation  can  be  obtained  accordingly.  Using  the 

pseudo  block  but  delete  the  row  and  column  where  K  situated  from 

D 

A(S),  then  the  determinant  of  one  less  order  can  be  obtained 

AmCS)*! A(S.  delete  the  row  and  colume  where  Kjj  situated! 

=MGtS)MH(S)“=0  (3.2.2) 

and  the  number  of  the  zero  points  and  their  quantity  of  the  open- 
loop  can  be  obtained.  From  A^US)  and  ,  the  transfer  function 

-’MS)  of  the  open-loop  can  be  obtained 

=  '3-2-3' 

based  on  open-loop  frequency  ^(jw-  and  stability  analysis,  the  polar 
coordinate  (Nyquist)  diagram,  logarithm  amplitude  diagram,  phase- 
frequency  (Bode)  diagram,  logarithm  amplitude-phase  (Nichols,* 
diagram  can  be  generated  and  the  s^'ability  and  width  of  the  POGO 
system  can  be  determined.  (See  figures  4.2,  4.3.) 

IV.  Single-Transfer  and  Estimation  Method  of  POGO  Stability  Analysis 
POGO  has  its  own  special  features.  In  fluid-solid  coupled 
dynamic  analysis,  some  systems  req*aire  at  least  two  modular 
vibrational  states  to  cause  the  unstable  vibration  of  the  entire 
system.  For  example,  the  gas  dynamic  vibration  of  aircraft  wings 
require. bending  and  twisting  vibrational  modes.  Compared  with  the 
inertia  and  weight  of  the  cross  section  of  the  structure,  if  the  gas 
dynamic  force  is  large  enough,  then  a  larger  deviation  of  the 


key :  1  -  pseudo  box 


2  -  effective  thrust 


3  -  post-pump  pipeline  and  thrust  chamber 

4  -  pump  and  longitudinal  oscillation 

5  -  compensation  pipeline 

6  -  pressure  retaining  pipeline 

7  -  compensation  pipeline 
S  -  five  pass 

9  -  main  line 

10  -  liquid  reservoir 


Fig  4.1 


.^0 


characteristic  frequency  of  the  system  can  be  caused.  However,  for 
some  other  systems  such  as  galloping  or  pogo.  the  fluid  or  pulse 
dynamic  force  is  relatively  small  compared  with  those  of  larger  and 
heavier  structures  and  only  small  deviation  of  the  system 
characteristic  frequency  can  be  caused.  Tlie  stability  analysis  of 
these  systems  can  be  carried  out  with  single-valence  modular  state 
parameter  method.  In  this  section,  we  consider  only  the  single¬ 
valence  modular  state  of  the  rocket  structure  and  the  propellent 
system  and  use  them  in  the  POGO  stability  analysis.  The  number  of 
equations  involved  in  greatly  reduced  in  this  way.  Based  on  the 
general  feature  of  the  structure  of  pressure  retainer  (see  figure 
4.1),  19  equations  for  the  19  parameters  can  be  obtained  from  the 
discussion  in  section  II.  The  stability  analysis  can  be  carried  out 
using  the  method  discussed  in  section  III;  using  either  the 
characteristic  root  of  the  close-loop  equation  or  the  transfer 
function  based  on  pseudo  block  principle.  One  can  also  carry  out 
parameter  elimination  from  the  .19  equations  and  the  transfer 
equation  of  the  loop  in  question  '.see  table  4.1)  is 


A'-(7(.S)r 

whereas  the  transfer  equation  of  the  feedback  loop  is 


(4.1) 


7’- //(A)  A 

Then  the  transfer  function  of  the  close- loop  with  output  and 
feedback  is 


F(S)  =  t?(6')/Cl-<7(5)H(S)] 


XI 


Table' 4.1:  Transfer  function  corresponding  to  figure  4.1b 


key:  1  -  tank  bottom  total  outlet  cross  sectional  area  /  N 


T I  "‘ch(tri)eh(sr,')+lttL,sh(STi)sh(sTi)2/lr,L,N>) 

Tt  ««Cr,cA(»r,)*A(*r,)]/L,+[r,jA(ir,)cA(4r,)]/(£,,/V») 

T  a  *cA(*r,)+[r,  t,cA(iri)*/i(jr,)]/[raLtAr'j/i{ir,)] 

Tt  -Cr.*A(ir,)]/£„+[r,cA(ir,)]/[LuV'+A(*r,)] 
r,,-[cA(*r,)+(r,L,tA(*r,))/r,]r,+C7’:L,jft{sr,)]/r, 

T  ,o~Lch(sr,)+(Y  ,Ltsh{su))/T^2T,HT,L,sh(.sTt)yr.. 

T  o^-[(r,tA(*r, ))/£,, +y,cA(*r,)]7',  +  rtcA(ir,) 

T  oo=[(TjtA(4r,))/£,,+y,cA(*r,)]7'i+r4cA(ir,) 

7'^'o’“[7’,ot,*A{4r,)]/r„  To'o^Croo  t,jA(4r,)]/r, 

iVi(*)-C5P.(r,,ro'o-r/oro^)]/^,(e) 

i>,(*)-r/o(SC,+(m+l)/2,)  +  ro'o 

Hi(*)-C(m+n2,/Z.-a<,^.(p)/(S,^.{e))]/VS,W,(0/I>.{i) 

•  Ht{s)^-sN  Pn’KM/DM 

■  H,(s)-m2sNA,'P.i.(t)/4.He),  Z,=Z,+Z4+Z<. 

HM~-N  AtC,N,(.s),/DM,  F(<)-Z,(m+»/Z, 

HtU)’^-C,H,(s)/A\ 

C,-Cma.^.(/p)^^'S*]/[a4,^,(p)F/.(*> 

C,-[m,^,»«p)^,.]/2pA,^,(0^,(/)F,(*)3 

P.m^9K'=*Th,^,{f)lg 

At’-A.'^-m&mamwN  ® 

AT'  -1*4 
Ff(j)=S'+2{AWf5+w^* 

H(i)  =H,(*)  +  Hi{*)  +  //,(4)  +  H.(i)  +  W,(j)  =  H,{4)  +  ;7/:(.0 
G(i)  »ff,(e)*/(i*+2{,ii),t+w,*) 

-C.(e)//,(«,)/(2«.) 


and  the  open-loop  transfer  function  is 


«s)=(?(s)«s).r^,,  t'-*!..  (4,4. 

where  X.  T  are  the  velocity  and  equivalent  thrust  of  the  engine  at 
its  normal  state.  The  expressions  for  <j{S)  and  H(S)  are  shown  in 
table  4.1. 

Based  on  Norquist  stability  criteria,  the  system  stability  can 
be  determined  through  the  calculation  of  the  characterist ic 
frequency  ^fjw)  of  the  open-loop.  From  table  4.1  one  has 


then 


^'(y<»){C<r,(e)cos0c3/2£.<».)}e»*c 

fic=/p-T(6)i-a>0/(2f^03 
//(/<»)=  I  i/(/(a)|e'*« 


^(;o)= 


G^e)\  Hiio)  I  cc60c^,t^ 

2^mC0it 


■  4 .  ?  i 


If  this  system  exhibit  iwj;  l  tor  :^{jw):.-li  for  any  phase 
orientation,  then  this  system  is  stable  (or  unstable). 

If  there  exists  and  i04(0,)-O, 


0e(at)- 


and  we  have 


(4.10) 


l^cl  =  =<?.(«)//»(«•.».)/(  2?^.) 


(4.11) 


(4.12) 

based  on  general  control  theory,  the  stable  amplitude  width  AdB) 
of  the  system  POGO  is 

A((/yn  =  20log  .,(l^(/a),)|)-'=  -201og  „a./^)-20log  ..(Wfc.),(Sr>l))  (4 

where  Ct.  has  the  same  meaning  as  before. 

If  the  open- loop  frequency  is  maximum  at  a=o„  and 

I0(;cu)l.«--=l^(;o)*)|  =  1«4»|>1  (4,14) 

then  and 

■  ■  ,4  15) 

and  we  have 


OaM  =  c»s-'C2^<b,/CC?„»  1  //(/«, )  I )] 


(4.16) 


then  the  phase  width  is 


■  6r  =  9iiia)r)  +  6c(.cOr) 


(4.17) 


If  system  |0cl<l  and  |<P«|.''l  at  non-sero  phase  orientation,  then 
this  is  called  the  phase-orientation  stable  system  and  pure  phase 
shift  would  cause  unstable  POGO  of  the  system. 

So  far.  except  for  the  assumption  of  single  modular  state  and 
single  propellent,  there  is  no  other  assumptions  made  for  the  POGO 
stability  analysis.  From  equations  (4.12)  and  (4.13)  we  knov?  that 
1.  the  rocket  body  structure  damping  ratio  A  is  small: 


2.  the  mass  is  somewhat  large  or  i.«,j  is  small: 

3.  rocket  body  gain  high;  and 

4.  the  real  part  H  ( cocj  of  feedback  transfer  function  is 

1 arge . 

All  the  factors  mentioned  above  will  decrease  A(dB)  and  unstable 
POGO  of  the  system  is  even  more  easy  to  occur.  This  method  and  its 
conclusions  also  apply  to  other  stability  analysis  of  similar  fluid- 
solid  coupling  systems. 

From  the  results  of  the  root  of  characteristic  equation  and  the 
above  analysis,  one  can  realise  that  the  frequency  for  POGO  to  occur 
is  very  similar  to  the  modular  state  frequency  of  rock  body 
structure.  In  other  words,  if  POGO  should  occur,  the  frequency 
should  be  around  the  structure  single-valence  modular  frequency.  One 
can  also  know  that  the  feedback  transfer  function  of  POGO  system  has 
its  own  characteristics:  namely,  before  the  measures  were  taken  to 
overcome  the  unstable  pogo.  the  characteristic  frequency  H(jw)  is 
slow-changing  with  respect  to  Gfjw)  and 

|/f(/0Je)|  |//(/w,)|  «  l/-/(;(y„)|  =;  |//(;w,)| 

then  we  should  obtain  the  following  estimation  equations 


A(c/S)  =  20  log  „(S./L) 


(4.20) 


key:  1  -  polar  coordinate  diagram 


2  -  logarithm  amplitude  -  phase  diagram 

3  -  logarithm  amplitude,  phase  -  frequency  diagram 

4  -  tc‘  ~  fi^equency  diagram 


Fig.  4.2: 


Single-transfer  stability  analysis  of  pressure 
retainer  free  system  {C_=0) 

ci 


key:  1  -  polar  coordinate  diagram 


2  -  logarithm  amplitude  -  phase  diagram 

3  -  logarithm  amplitude,  phase  -  frequency  diagram 

4  -  Jc*  ®r  “  frequency  diagram 


IfS.S  11/  lil.l  l».6  »l  .!U.: 


Fig.  4.3:  Single-transfer  stability  analysis  of  system  vnth 
pressure  retainer  ( Cq=18 . 57cm2/kg) 


fable  4.2 r Assumed,  pafametefs  for  space  vehicle 
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key:  1  -  fixed  parameters 

2  “  system 

3  -  oxidizing  agent 

4  -  density 

5  -  resistance  of  thrust  chamber 

•  •Ir.i'd---  cross  sectipnal  area  of  suction-  pipe..^ 


7  -  length  of  suction  pipe 


8  •*  speed  of  sound 


9  -  resistance 

10  -  inertia 

11  -  modular  ratio 

12  -  thrust  chamber  time-delay  constant 

13  -  modular  damping  ratio 

14  -  (b)  time  dependent  parameters 

15  -  system 

16  -  oxidizing  agent 

17  -  percentage  of  time  of  flight 

18  -  modular  frequency 

19  -  macroscopic  mass 

20  -  pump  dynamic  gain 

21  -  liquid  height  of  liquid  reservoir 

22  -  pump  gas  softness 


.  ♦  ' 


.  ^8 


Table  4.3 


key:  1  -  related  parameters  and  pressure  retainer  parameter 
in  table  4.2  for  95%  time  of  flight 

2  -  POGO  stability  equation  and  calculated  results 

3  -  estimation  method 

4  -  single-transfer  method 

5  -  pressure-retainer  free 

6  -  pressure-retainer 

7  -  note:  Effect  of  tank-bottom  flov;  is  not  included. 
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where  S((/5),  are  the  amplitude  of  the  POGO  system  and  the  estimated 
value  of  the  phase-stability  width.  Following  this  method,  the 
varying  scanning  frequency  a  can  be  manually  fixed  at  the  single¬ 
valence  modular  frequency  o,,.  If  this  algorithm  is  followed,  than 
the  amount  of  calculation  is  reduced  and  an  overall  estimation  can 
be  accomplished. 

Literatures  [4]  and  [5]  have  studied  the  estimation  method, 
however,  the  error  and  direction  of  deviation  were  not  given  and 
there  these  methods  can  not  be  used  independently,  in  table  4.2.  a 
set  of  aircraft  parameters  were  listed  (see  [5])  and  table  4.3  is 
the  calculational  results  (with  and  without  pressure  retainer)  based 
on  parameters  in  table  4.2.  In  these  calculations,  single-transfer 
and  estimation  methods  were  used.  Figure  4.2  and  4.3  are  the  polar 
coordinate  (Nyquist),  logarithm  amplitude-phase  (Nichols),  logarithm 
amplitude,  phase-frequency  (Bode),  and  S,.(f)  and  ( f) -frequency 
diagrams  corresponding  to  the  single-transfer  method.  From  these 
diagrams,  the  width  of  the  system  stability  and  the  error  of  single¬ 
transfer  and  estimation  methods  can  be  deduced. 

From  these  tables,  figures,  and  calculation  experience,  it  is 
known  that  before  measures  were  taken  to  overcome  POGO  of  space 
vehicles  (such  as  installation  of  pressure  retainer),  estimation 
method  provides  good  accuracy  (relative  error  smaller  than  3dB) . 
however,  after  overcome  measures  are  taken  and  estimation  method  is 
used  to  select  parameters,  the  error  is  large  (a  few  tens  of  dB)  and 
the  results  tend  to  be  not  as  reliable.  Therefore,  more  accurate 
methoas  such  as  single-transfer  or  the  methods  explained  in  section 
III  should  be  used  to  check  the  results. 


Another  preliminary  estimation  method  for  the  analysis  of 

stability  is  to  use  the  characteristic  of  open-loop  frequency 

4>(  jw)=G(  jw)H(.iw)  .  When  at  zero  phase  orientation,  if  w-w  than 

c 


(4.22) 


then  the  POGO  system  is  not  stable.  If  is  in  the  vicinity  of  the 
peak  frequency  w*  of  !H(jw)l,  stability  of  the  system  can  be 
determined  from  the  closeness  of  w*  and  Wj^.  Therefore,  the  first  of 
determining  the  stability  of  the  system  is  to  obtain  the 
characteristic  frequency  w*  of  H(jw).  To  make  the  problem  simpler, 
one  can  consider  the  five  cross  sections  t-l-3-5-D  of  figure  4,1, 
there  are  four  transfer  functions  (from  section  II) 
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Z,  ) 

based  on  the  "open"  (Pi=Q,  or  "close"  Q^=0)  at  cross 

sections  t  and  D,  one  can  obtain  the  low-frequency  and  high- 
frequency  resonance  equation  and  formulae  for  various  feedback 
systems  (mainly  the  pipeline  and  pump  systems). 


V,  Analysis  and  Discussion  of  Major  Parameters 

The  analysis  in  this  section  is  based  on  the  related  equations 
in  section  IV  (see  table  4.1). 

1.  outlet  flow  from  the  bottom  of  the  tank 


3/ 


Table  4.4 
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1  -  component  parameters 

2  -  main  component  of  the  system 

3  -  long  suction  pipe 

4  -  long  suction  pipe,  pump 

5  -  long  suction  pipe,  pressure  retainer,  pump 

6  -  long  suction  pipe,  pressure  retainer,  short  suction 

pipe,  pump 

7  -  boundary  condition 

8  -  open-close 

9  -  open-open 

10  -  open-close 


13  -  open-opeu 

14  -  open-close 

15  -  open-close 

16  -  open-open 

17  -  close-close 

18  -  open-close 

19  -  open-close 

20  -  open-open 

21  -  close-close 

22  -  pump 

23  -  softness 


11  -  open-open 

12  -  open-close 


Table  4.4  (key  continued) 

24  -  resistance 

25  -  inertia 

26  -  short  tube 

27  -  inertia 

28  -  pressure  retainer 

29  -  softness 

30  -  inertia 

31  -  long  pipe 

32  -  speed  of  sound 

33  -  See  figure  4.1  and  section  II  for  related  transfet 

function.  If  suitable  boundary  conditions  are 
plugged  in.  the  following  frequency  equations  can 
be  obtained. 

34  -  gas  bubble  frequency 

35  “  high  frequency  condition 

36  -  high  frequency  condition 

37  -  reverse  resonance  fi'eauency 

38  -  Long  pipe  inertia  i.e  Li  and  lencrth  is  li. 


When  calculate  or  test  the  modular  states,  the  tank  bottom  is 


either  closed  or  the  liquid  is  not  flowing.  When  the  flight  POGO 
analysis  is  carried  out,  the  tank  bottom  is  opened  and  the  liquid  is 
flowing  and,  therefore,  the  contribution  of  the  outlet  flow  from  the 
bottom  of  the  tank  to  the  macroscopic  force  should  be  considered. 
This  is  reflected  in  the  H2(S)  and  H3(S)  in  table  4.1.  This  factor 
will  contribute  to  the  stability  of  the  system  and  it  is  over¬ 
conservative  if  this  factor  is  not  considered. 

2.  local  pump  longitudinal  \  ..bration 

The  effect  of  local  pump  longitudinal  vibration  is  reflected  in 
the  contribution  to  the  macroscopic  force  and  is  manifested  in  Ho.CSl 
and  H4(S)  terms.  When  the  frequency  of  the  local  pump  longitudinal 
vibration  Wp  is  close  to  or  overlapping  with  the  rocket  body 
structural  modular  state  frequency  Wj^,  installation  of  pressure 
retainer  will  not  suppress  POGO.  As  a  result,  it  is  very  important 
to  clarify  the  parameters  for  local  pump  longitudinal  vibration 
before  flight. 

3.  composite  pre-ump  pipeline  resonance  frequency  w* 

Experience  shov/ed  that  when  the  characteristic  pre-purnp 

frequency  w*  is  close  to  the  rocket  body  structure  modular  frequency 
w^.  unstable  POGO  of  the  system  will  occur.  Hence,  to  separate  w* 
far  from  Wj^  is  one  of  the  important  criteria  for  the  design  of  space 
vehicles.  (9)  This  can  be  made  clear  from  the  following  analysis. 

i.  When  the  tank  bottom  outlet  flow  and  local  pump  longitudinal 
vibration  are  not  considered,  we  have 

HiS)  =  H,iS)=\EiS)-^]NS,A,t,SK.iS),V(,S)/rKS) 

t  ( 5 .  J .  i ) 
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the  meaning  of  these  symbols  is  explained  in  table  4.1  and  figure 
4.1  (same  for  the  following  equations). 

ii.  When  there  is  only  one  pipeline  section  and  pres-sure 
retainer  prior  to  the  pump,  the  E>(S)  in  the  above  equation  takes  the 
form  of  D(iw)  and 


Ltg(;r6)/(2o)i))  l+;C.«/?uf<)  -Cjaco  /J 

+  /(m  +  l)/,o) 


(5.3.2) 


where  is  the  composite  softness  of  the  gas  bubbles  in  the  liquid 
within  the  pipeline. 

iii.  Wlien  the  inertia  of  the  pressure  retainer  and  resistance 
are  not  considered,  we  have 


D(y<»)=Z,(j<a) 


go/(2<a,) 


Where car«a,/(2/i>and  is  the  first  valence 
pipeline,  and  «&  =  C/i(C,+Cp+C,)]-''- 

iv.  Let 


(5.3.3) • 

frequency  of  the  open-close 
is  the  gas  bubble  frequency. 


p/  ;r<i)/(2c?i) _ a’ 

tg(W(2o).))  col  (5.3.4) 

if  w-w* .  then  F(w*)=0  and  is  called  the  composite  resonance 
frequency  of  the  pipeline.  When  D(jw*j  is  a  minimum.  H(jw*)  is  a 
maximum  and  S.fw*)  is  also  a  maximum.  Therefore,  when  w*»Wn, 
unstable  POGO  of  the  system  would  occur  very  easily.  However,  it  is 
not  true  that  POGO  should  alv/ays  occur  because  there  are  other 
parameters  which  also  control  the  mechanism  (see  below) . 


4.  Overlapping  of  w*=Wn 


■From  above,  when  1)  the  tank  bottom  outlet  flow  is  not 
considered.  2)  local  longitudinal  pump  vibration  occurs.  3)  no 
pressure  retainer,  and  4)  only  one  section  of  equivalent  pipeline 
exists  prior  to  the  pump,  we  can  have 


=<?,(e)/f  «(©.)/)  2o),''= 


G»ie)NS.A%<l>0Rj' .  ,  A»^/r,{a, _ f 


From  the  above  equation,  one  can  learn  that  even  if  w*“Wn.  is 

possible  under  certain  conditions  and  pogo  will  not  occur.  Among 
these  conditions,  the  gain  of  the  engine  Eq  is  an  important 
parameter.  The  significance  of  eq.  (5.4.1)  is  that  the  effect  of 
related  parameters,  modular  states,  or  structural  gains  on  the 
system  stability  can  be  determined. 

5.  parameters  and  location  of  pressure  retainer 

So  far  there  are  three  methods  to  select  the  parameters  for 
pressure  retainer 

i.  reverse  resonance  frequency  method,  namely.  w^^=*w*-=Wj^: 

ii.  POGO  stability  dependability  method:  and 

iii.  compromised  method,  vzhich  suppresses  the  unstable  POGO  of 
the  system  and  at  the  same  time  keep  the  typical  frequency  of 
pressure  retainer  line  to  different  from  the  fixed  frequency  w^w^. 

From  table  4.4.  the  effect  of  the  parameters  and  location  of 
pressure  retainer  on  the  w*  can  be  observed.  If  the  pressure 
retainer  is  placed  near  the  pump  inlet,  the  efficiency  of  changing 
w*  is  enhanced  and  the  adjustment  of  pressure  retainer  parameters  is 
made  easier,  while  on  the  other  hand,  the  water  impact  pressure  is 
reduced  and  the  pressure-reduction  due  to  pre-pump  high  pulse 
pressure . 
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From  table  4.4.  it  is  known  chat  after  the  pressure  retainer  is 
added,  when  the  pre-pump  pipeline  is  under  the  condition  of  open- 

A  A 

close,  the  first  and  second  valence  frequencies  w  ^  and  w  are 


o»r-l/ v//.(C:+C,)  (5.5.1) 

(^^i=y/c~+c;/y/u;+n)c7r, 

(5.5.2) 

Attention  should  be  paid  to  the  range  of  applicability  of  above 
equations.  Generally.  can  be  reduced  by  increasing  C^.  however. 
w*2  not  be  increased  infinitely  by  reducing  (Ig+la)  the  w*2 
will  be  lower  than  the  original  second  valence  frequency  (without 
Sddition  of  pressure  retainer) . 

6.  equivalent  speed-of-sound  of  pipeline  liquid 
From  (5.3.3),  it  can  be  learned  that  a^  first  affect  w,  then 
w*,  and  finally  the  P060  stability. 

Not  only  is  ae  related  with  the  property  of  the  liquid,  it  is 
also  related  with  the  gas  content  Vg  generated  during  the  flow 
process,  pipeline  material,  and  boundary  supporting  conditions.  When 
the  volume  of  liquid  in  the  pipeline  is  Vj.  the  total  volume  is 
V^VjVg  and  the  composite  density  p  is 


P^PiVilV^P,V,lV 

and  the  composite  elastic  modulus  K  and  a^  are 


(5.6.1) 


=  A',/[l  +(K,/k)(A:,/A:,  -1.)]  (5.6.2) 

a.=  v/A/p"  ly/i^(KIE)CDl6TC 

(5.6.3) 

where  D.  5,  E.  p  are  the  inner  diameter,  thickness,  elastic 


modulus,  and  Poisson's  ratio  of  the  pipeline  material.  Note  that  C 
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is  related  with  the  boundary  supporting  condition.  When  the  upstream 
end  of  the  pipeline  is  fixed  and  the  downstream  end  is  free.  C-1- 

/i/2;  when  there  is  no  axial  strain  of  the  entire  pipeline,  C=l-/x’i 
when  there  is  not  axial  stress.  C»l. 

The  gas  content  in  the  pipeline  liquid  Vg/V  and  the  theoretical 
and  experimentally  obtained  a^  curves  are  shown  in  figure  5.1.  tlO] 

7.  pre-pump  liquid  gas  content  softness  and  pump  gas 
softness  Cp 

From  (5.3.3),  the  major  factors  influencing  w^  is  and  Cp. 
Theoretical  determination  of  these  quantities  is  difficult  and  only 
experimental  method  can  be  applied.  Based  on  the  experimentally 
determined  w*  and  and  through  equation  (5.3.4)  (F(w*)-0).  wj^  can 
be  obtained  and  Ci+Cp  can  be  calculated. 

8.  boundary  condition  of  pre-pump  pipeline 

Generally,  when  the  characteristic  frequency  of  the  pre-pump 
pi,.'eline  is  calculated  or  experimentally  determined,  the  effect  of 
the  boundary  supporting  conditions  (both  ends)  is  significant. 
However,  in  POGO  calculation,  the  boundary  condition  of  the 
connection  between  pipelines  is  not  clearly  defined.  This  is  because 
the  boundary  condition  should  be  determined  by  th-  continuous 
condition  of  the  system  loop.  For  example,  when  Ca+Cp+Ci»0,  from  eq. 
(5.3.4)  F(w*)-0 

ca*=ffo,(2n-l)/(2/,),  n=U  2-  (5  8  1) 

and  this  corresponds  to  the  open-close  boundary  condition. 

Similarly,  when  Ca+Cp+Ci  «o6'.  then 


0*~nnajltf  n=®l,  2,  ••• 

and  this  corresponds  to  the  open-open  boundary  condition.  In 
reality. 


0<C.+C,+C|<co 


(5.8.3) 


and  if  F(w*)“0  is  used  to  calculate  composite  resonance  frequency  of 
the  pipeline,  then  we  should  come  up  with  the  boundary  condition  for 
either  open-open  nor  open-close. 

9.  dynamic  gain  (m+1)  of  the  pump 

Dynamic  gain  of  the  pump  is  one  of  the  quantities  that  is 
somewhat  difficult  to  determine.  It  affects  the  gain  of  the  engine 
Eg  and  other  quantities  (see  eq.  5.4.1)  and  furthermore,  the  system 
stability.  One  of  the  key  parameters  in  POGO  calculation  is  m+l*. 

From  the  relationship  of  the  variables  in  section  II.  we  have 


(m  + 1)= /^,(;ta)/’o/CF»(2e(/o)  +  Zo(;ta)] 

(5.9.1) 

where  Zq.  Zq.  Zq  (see  table  4.1)  can  be  obtained  from  calculation. 
Pu  and  Pj)  are  the  pressure  pulse  before  and  after  the  pump,  v/hich 
can  be  obtained  from  experiment. 

10.  gain  of  the  engine  Eg 

From  eq.  (5.4.1).  when  Eg<A50p/ (SgOg) .  we  have  5*<0  and  the 
design  of  the  engine  is  adequate.  No  POGO  will  occur  even  if  the 
resonance  frequency  w*  of  the  pre-pump  pipeline  is  overlapped  with 
w^.  This  is  also  one  of  the  most  important  design  criteria  for  the 
engine  of  space  vehicles. 

11.  gain  of  rocket  body  structure 


^0 


The  structure  gain  G^(e)  can  be  defined  based  on  eq.  (2.1.3) 
and  from  eq.s  (4.12)  and  (5.4.1),  it  can  be  seen  that  or  its 
estimated  value  L  are  proportional  to  •  Since  Gjj(e)  can  be 

expressed  in  terms  of  the  time-of-f light  as  Gn(e,t).  therefore,  the 
trend  of  ^e'(t)  and  ^c(t)  is  similar  to  that  of  6^(6. t)  and  the 
accuracy  of  the  result  can  be  determined.  Gain  of  rocket  body 
structure  is  also  one  of  the  key  parameters  in  POGO  design. 

Sometimes  even  slight  change  in  the  structure  modular  frequency  will 
result  in  significant  change  in  structure  gain. 

12.  structure  modular  damping  ratio 
.f,.  consists  of  three  parts:  fluid  damping,  material  damping, 
and  structure  damping.  It  can  be  learned  from  eq.s  (3.1.4),  (4.13) 
and  (4.20)  that  the  magnitude  of  Cn  has  a  direct  effect  on  the 
accuracy  of  POGO  stability  analysis.  Hence,  this  parameter  is  an 
important  parameter  and  is  generally  provided  by  test  or  experience. 
The  test  result,  even  if  identical  method  were  used,  will  show 
discrepancy.  The  discrepancy  of  the  f,.  of  same  order  based  on 
various  methods  will  be  even  greater.  In  reference  (12],  was 
obtained  based  on  various  test  and  processing  techniques  and  the 
discrepancy  is  very  significant.  From  a  broader  prospective, 
decreases  with  modular  frequency  as  shown  in  figure  5.2.  [12]  If 
high  vacuum,  low  gravity,  and  other  conditions  that  can  not  easily 
simulated  on  earth  is  considered,  it  is  even  more  obvious  that 
dependability  analysis  of  Cn  should  be  carried  out.  Because  of  the 
difference  in  test  methods,  boundary  supporting  conditions,  payload 
conditions,  and  the  non-linear  effects,  the  obtained  from  the 
full  scale  test  will  show  large  discrepancy  and  is  not  dependable. 

V/ 


key:  1  -  sinusoidal  fixed  frequency  test  data 

2- 4  excitors  random  test  data 

3- 3  excitors  random  test  data 

4  -  90%  of  the  data  are  lower  than  this  line 

5  -  modular  damping  ratio 
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Fig.  5.2:  Modular  damping  ratio  obtained  from  various 


methods 


Even  though  in  the  stability  calculation,  the  value  of  is  taken 
to  be  0.01  and  0.007  by  groups  in  USA  and  France,  respectively,  it 
is  generally  believed  that  it  is  not  over-conservative  even  af  a 
lower  value  of  is  used.  If  the  suppression  of  POGO  is  difficult 
using  other  methods,  the  increase  in  damping  ratio  of  the  key  parts 
of  the  rocket  structure,  the  modular  damping  overall 

rocket  structure  will  be  increased  and  POGO  can  be  suppressed.  [13] 

13.  comparison  of  the  importance  of  various  structural  modular 
parameters  of  different  orders 

The  importance  of  the  structural  modular  parameters  of  certain 
order  is  evidently  determined  by  the  actual  problems  encountered. 

For  example,  if  the  structural  response  to  the  acute  bottom  seat 
vibration  (for  aircraft  or  satellites.'  is  considered,  the  modular- 
state  related  parameters  such  as  "effective  mass”  is  an  important 
basis  for  the  determination  of  the  importance  of  certain  modular 
states.  For  the  POGO  problem,  the  modular  states  which  cause  POGO 
are  obviously  the  important  states.  They  are: 

i.  large  mass,  low  frequency  Wj.^: 

ii.  higher  structural  gain  G^(e): 

iii.  small  modular  damping  ratio  and 

iv.  modular  states  which  make  w*  close  to  Wp. 

VI.  Dependability  of  POGO  Stability  Calculation 

The  design  of  space  vehicles  is  based  on  the  nominal  values  of 
various  parameters.  In  reality,  all  the  parameters  will  deviate  from 
the  nominal  value.  The  test  value  of  a  certain  parameter  will  also 
depend  on  the  location,  time,  installation  method,  temperature. 


humidity,  and  status  of  the  instrument.  Furthermore,  the  result  will 
be  different  if  test  methods  is  different. 

Since  the  development  cost  of  space  vehicle  is  huge,  to 
guarantee  flawless  design,  the  dependability  calculation  and 
analysis  is  important.  If  a  new  product  is  involved,  the  test  data 
may  be  scarce,  the  correct  data  may  be  difficult  to  obtain  from 
theoretical  calculation  or  from  a  single  test.  Often  it  would 
require  a  lot  of  man-power,  a  long  preparation  period,  and  a 
tremendous  budget  for  a  particular  test  which  may  be  extremely 
contaminating  to  the  environment.  For  the  systems  which  include 
parameters  of  this  kind,  it  is  better  to  use  the  computer  to  carry 
out  a  Monte  Carlo  system  dependability  analysis. 

From  sections  II  and  III.  it  is  clear  that  there  are  many 
primitive  parameters  involved  in  the  dependability  analysis  and  each 
parameter  has  its  own  probability  distribution.  ITie  combination  of 
all  the  parameters  represents  a  simulated  flight  and  the 
dependability  analysis  for  this  particular  flight  can  be  carried 
out.  Hence,  the  combination  of  random  sampling  of  all  the  parameters 
(based  on  its  probability  distribution)  and  the  estimation  criterion 
for  unstable  POGO  will  generate  n  sets  of  parameters  which  represent 
n  times  of  simulated  flight.  Tlie  POGO  stability  analysis  for  these  n 
flights  can  be  proceeded  and  than  the  probability  of  unstable  POGO 
can  be  calculated  based  on  some  credibility  criterion. 

There  are  constants,  time-dependent  constants,  random  constants 
of  various  probability  distribution,  and  time-dependent  random 
constants  in  a  POGO  system.  For  example,  modular  state  damping  ratio 
f,.  exhibits  the  minimum-value  distribution,  pump  dynamic  qain  (mi-1) 


exhibits  the  maximum-value  distribution,  pump  gas  softness  Cp  and 
pipeline  equivalent  speed-of-sound  a^  exhibit  normal  distribution. 
The  sampling  method  for  these  parameters  are: 

f./=i/.+5,InInCl/(l-i?.i)],  /=1,  2,  -n  fg 

(m/+l)=C/*-5*InlnCl/i?»/3»  ;  =  1»  2,  — »  (6.2) 

;=1,  2,  •••«  (6,3) 

where  U^.  U]<-  ^n-  determined  by  the  estimated  values  of 

related  overall  parameter  and  Rm  and  Rpj  are  the  number  of  uniform- 
divisions.  fi-Q  and  ffo  are  the  averages  of  related  overall  parameters 
and  estimated  value  of  standard  deviation.  S^j  is  the  normal  state 
random  number  and  j  is  the  designation  of  the  random  state. 

From  section  V.  one  learns  that  some  of  the  parameters  are 
important  but  difficult  to  .determine .  If  sampling  is  carried  out  for 
all  the  parameters,  and  then  POGO  stability  analysis  conducted  for 
all  the  combinations  (while  the  flight  condition  is  changed  in  unit 
of  second),  the  amount  of  calculation  will  be  significant.  From 
experience,  some  of  the  undetermined  parameters  can  be  manually 
reduced  or  combined  to  reduce  the  amount  of  calculation.  Likewise, 
the  number  of  combinations  can  be  reduced  or  the  flight  simulation 
can  be  carried  out  based  on  a  longer  time  unit. 

VII.'  Criteria  for  POGO  Stability  and  Dependability 

It  can  be  seen  from  sections  III  and  IV  that  POGO  stability  is 
determined  by  the  amplitude  width  A(dB)  and  phase  width  0'^. (degree)  . 

The  criterion  for  POGO  stability  is  similar  to  the  control 
system;  namely,  amplitude  width  A>6dE  and  phase  width  !  0^:  >  30'^'. 


This  corresponds  to  the  structural  safety  consideration  of  design 
safety  coefficient  of  2.  However,  this  stability  width  can  not  be 
extended  to  uncertain  parameters  such  as  Cp.  or  (m+i)  and 
dependability  statics  should  be  employed.  Under  these  circumstances, 
the  criterion  for  POGO  stability  is  3^:  namely,  the  probability  of 
unstable  POGO  as  a  result  of  the  random  combination  of  all  the 
parameters  is  less  than  0.00135. 

From  above,  it  can  be  concluded  that  a  better  dependability 
criterion  of  POGO  stability  is  for  the  nominal  widths  to  be  A.>6dB 
and  1  >  30°  in  addition  to  the  acceptable  probability  of 

unstable  POGO  of  P< 0.00135. 

VIII.  Conclusions 

1.  When  the  space  vehicle  engine  pulse  thrust  frequency  and  the 
rocket  body  structural  frequency  (w*  and  Wp  or  w^^j  are  near  or 
overlapped  and  the  structural  gain  is  sufficient,  then  the  pulse 
motion  of  the  liquid  in  the  fuel  tank  is  equivalent  to  the  pulse 
excitor  and  it  becomes  as  a  source  of  excitation  for  the  thrust 
system  behind  the  tank  bottom.  This  constitutes  a  positive  feedback 
system.  The  simplest  model  for  this  system  is  the  free-free  double 
mass-spring-damper  system.  This  model  connects  the  POGO  stability 
damper  analysis  principles  based  on  various  methods. 

2.  After  the  dynamic  equations  of  POGO  system  are  introduced, 
the  POGO  matrix  analytical  method,  computer  calculation  of  open-loop 
transfer  function.  POGO  single-transfer  method  and  estimation  method 
are  explained.  These  methods  can  be  used  in  various  design  stages 
and  can  be  complementary  to  one  another.  Through  the  use  of 


examples,  accuracy  of  the  estimation  method  is  evaluated  before 
suppressing  measures  were  taken.  However,  after  suppression  measure 
parameters  are  selected  for  during  selection),  the  application  of 
estimation  method  should  be  cautioned  and  calibration  with  more 
accurate  methods  should  be  performed. 

3.  Analysis  and  discussion  of  some  important  parameters  is 
given  in  section  V.  This  section  also  points  out  some  of  the 
important  parameters  in  the  design  of  space  vehicles,  engine  system, 
or  the  suppression  measures  of  unstable  POGO.  Furthermore,  the 
various  measures  for  suppressing  POGO  is  discussed. 

4.  Dependability  analysis  should  be  carried  out  for  POGO 
stability  analysis.  Based  on  the  test  data  and  test  experience,  the 
computer  simulation  (such  as  Monte  Carlo  method)  can  be  used  for  the 
dependability  analysis.  Tne  criterion  for  the  POGO  stability/ 
dependability  is  such  that  the  nominal  width  A;>6dB  and  : 

plus  the  probability  of  unstable  POGO  P< 0.00135. 
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